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The temperature dependence of the interlayer electrical and thermal resistivity in a layered metal are calcu- 
lated for Fermi liquid quasiparticles which are scattered inelastically by two-dimensional antiferromagnetic spin 
fluctuations. Both resistivities have a linear temperature dependence over a broad temperature range. Extrapola- 
tions to zero temperature made from this linear- T range give values that appear to violate the Wiedemann-Franz 
law. However, below a low-temperature scale, which becomes small close to the critical point, a recovery of 
this law occurs. Our results describe recent measurements on CeCoIns near a magnetic field-induced quantum 
phase transition. Hence, the experiments do not necessarily imply a non-Fermi liquid ground state. 
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Strongly correlated electron materials exhibit a subtle 
competition between a range of ground states including 
metallic, insulating, superconducting, antiferromagnetic, and 
paramagnetic.l 1] Often the metallic states are distinctly differ- 
ent from the Fermi liquid state characteristic of elemental met- 
als. Heavy fermion metals are particularly interesting because 
they undergo quantum phase transitions, and have non-Fermi 
liquid properties near the quantum critical point. lit] For ex- 
ample, the material family CeMIns (where M=Co,Rh,Ir) can 
be tuned through transitions associated with antiferromagnetic 
or superconducting order by varying magnetic field, pressure, 
or chemical composition. Understandingthese systems has 
motivated significant theoretical eff'ort.lJ.lslla] In CeCoInj at 
ambient pressure a quantum phase transition between super- 
conducting and metallic states occurs as the magnetic field H 
is tuned through a critical value i/e-lii8»5J Recent transport 
datalfioll, from the vicinity of this critical point, display an 
extraordinary violation of the fundamental Wiedemann-Franz 
(WF) law of metals, and have been interpreted ! 11] as a possi- 
ble signature of a non-Fermi liquid ground state of CeCoIns . 

In this Letter we consider the WF law near a quantum- 
critical point with a goal to understand what observable WF- 
violation (WFv) reveals about the electronic ground state. The 
law states lll2ll that the electrical resistivity p{T) is equal to the 
electronic thermal resistivityfis'] w{T). While it should not 
hold at finite T, since inelastic scattering may be important, 
the WF law must be obeyed by a Fermi liquid at T = where 
scattering is due to static defects. 1 14] An intriguing aspect of 
the CeCoIn5 data is that T — > intercepts of p{T) and w(T), 
extrapolated from the lowest observed T of roughly 50 mK, 
have a finite diff'erence w{T ^ 0) - p{T — > 0) that increases 
as H is tuned towards He- This makes it appear that WFv 
might persist down to T = 0, revealing a breakdown of the 
Fermi liquid ground state. However, we show that a Fermi 
liquid subject to scattering by 2D critical spin fluctuations ex- 
hibits WFv at finite T, and in T — > extrapolations made from 
above a low-temperature scale, while still obeying the WF law 
at r = 0. Since a Fermi liquid model captures the distinctive 
H and T dependence of the CeCoIns data, these data[10] do 
not necessarily imply a failure of the Fermi liquid picture. 



We first place WFv within the context of transport phe- 
nomenology in CeCoIns. The resistivity p{T) is linear in T 
below 40 K over a wide range of field H and chemical doping, 
which could be an indication of quantum critical behavior |.15ll 
although its precise origin is unclear. ||3 For current along 
the stacked Celns planes, i.e., intralayer current, linear-T p(T) 
extends down to 5 K, roughly the same T below which anti- 
ferromagnetic correlations appear lll7[|l8|] ). before a downturn 
with decreasing T . Below 1 K, intralayer p{T) and w(T) con- 
verge, suggesting that the WF law is obeyed for intralayer cur- 
rents at sufficiently low T (the WF law is also approximately 
obeyed above 5 K in intralayer datafl^). Interlayer p{T) is 
T-linear down to the lowest measurable T with no trace of a 
downturn. For H well above p{T) and w{T) extrapolate 
to similar values at T =0. But as H is decreased towards 
He, the interlayer thermal resistivity w{T) undergoes a rigid 
upward shift that results in extrapolated T — > WFv. Since 
w{T ^ 0) > p{T — > 0), the WFv cannot be due to heat trans- 
port by neutral carriers Ii 19] but suggests instead that inelastic 
scattering contributes to the T — » resistivities. 

Since the T — > WFv is less robust than T-linear pij), it 
can be plausibly attributed to a diff'erent mechanism. WFv 
might result from scattering by inelastic antiferromagnetic 
spin fluctuations active at low temperatures, with T-linear 
p{T) determined by another scattering process important over 
a wider temperature range 1 6] . Taking into account its ex- 
perimental Fermi surface I20[ we model CeCoIns as a 
quasi-2D metal subject to strong two-dimensional antifeiTO- 
magnetic spin fluctuations and study interlayer transport us- 
ing a Boltzmann-equation description. The approach aims to 
understand the low-temperature WFv in interlayer transport 
without addressing other unusual properties associated with 
larger temperature scales. 

The relaxation rate of interlayer currents t"' = Tq' + f"' 
is written as a sum of an isotropic, elastic part t„' and an 



anisotropic, inelastic part f 



, e, r), coming from 



critical spin fluctuations. At low-T, in the critical region, 
will be active only near hot spots, i.e. near a pair of points 
on the 2D Fermi surface (of a single layer) that are connected 
by a spin-ordering wavevector Q. (To have finite interlayer 
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current, the Fermi surface must have some modulation along 
the interlayer momentum fe,. If spin fluctuations are spatially 
confined to a layer, and thus capable of imparting an arbitrary 
momentum transfer along L, then hot 'spots' will actually be 
lines nearly-parallel to the k,-axis on the quasi-2D Fermi sur- 
face.) The inelastic scattering rate depends on the direction (f> 
of electron momentum within a layer, on energy e and tem- 
perature T. The electrical and thermal resistivity can both be 
written, using rj^^iT) = p'^(T) and 7/2(7") = w"'(r), as 

,-,r, = ,;.£...(^)£;;|A«,,,n (1, 

where A-\(P,x,T) = 1 + TQf-\^,x,T), t]„^ = p-^\3 In^fl^ , 
po the zero-temperature interlayer resistivity and foix) = (1 -H 
e*)"' is the Fermi function and vertex corrections have been 
omitted (see Ref. ,23). Here, we assume symmetry-equivalent 
hot spots spaced by angle 4> with one at = (this assumption 
is not crucial, and any other distribution of distinct hot spots 
gives similar results). At T = we have t"' = Tq ' and p(0) - 
w(0) - Po, so the WF law is obeyed at sufliciently low T 
in this Fermi liquid model. But since our understanding of 
r = properties is based on measurements made at finite 
temperature, the effect of f ' on T ^ extrapolations of the 
model should be considered. 

The rate of scattering of electrons by spin fluctuations OlSll 
(the lowest-order electron self-energy with the spin suscepti- 
bility as the boson propagator) is 

T ((p,x,T)^2g, ) — — Xk-k'(kBT[x - X ]) (2) 

MX) 

where the k' sum is done in the usual way as an integral over 
linearized band energy x' = eu' /kgT and position on the Fermi 
surface with k - k' dependent only on <p,(p'. no(x) - 

(e-* - 1)"' is a Bose function and y'^ioj) is the imaginary part 
of the spin susceptibility, [4^ 2Ac i25J 

ikBT2)-\-^\a>) = -i-^ + + {kf^-^, (3) 

the energy scale ^^7^2 is the width of the spin-fluctuation spec- 
tral function at typical q (it is proportional to the parame- 
ter Tq of Ref. ^ and gs a coupling constant . The factor 
WqA:^. = (q ± Q)^ for |q + Q| « kf, while for |q ± Q| ^ kj, it 
is roughly q-independent. The spin fluctuations are assumed 
two-dimensional (the 3D case is discussed below) so Xni'^) 
is independent of q,. Also, we take for the spin-correlation 
length(25:] at T « T2 

i^qoT^ = r + CT/T2 (4) 

where c is a constant of order unity and r, the tuning param- 
eter, measures the proximity to the quantum critical point . 
For the magnetic field-tuned quantum critical point of inter- 
est r depends on H and vanishes at H = He- In the expres- 
sions above we have ignored logarithmic corrections, associ- 
ated with the system being in its upper critical dimension ll25ll . 



We discuss three temperature regimes for the spin fluc- 
tuations (always remaining close to the critical point where 
\r\ « 1), which are indicated respectively as I, II and III in 
the upper-left inset of Fig. [1] At T << T2 only spin fluctua- 
tions close to an antiferromagnetic Q are thermally excited so 
only electrons near hot spots encounter inelastic scattering. In 
the low-temperature region (I), defined by T « rTo « T2 
(with r > 0), the correlation length ^ is determined by the 
tuning parameter, r in Eq. |4] with temperature giving a weak 
correction. This may be distinguished from an intermediate 
temperature regime (II), rT2 « T « T2, in which this situ- 
ation is reversed. At high temperatures (III), T > T2, all spin 
fluctuations are thermally excited so there are no hot spots and 
f ' is independent of 0. 

Another significant temperature scale is that at which the 
inelastic scattering rate becomes comparable to Tq' (this is 
clarified below in the discussion of orbital effects of the mag- 
netic field). According to Eq. |2] the strength of inelas- 
tic scattering is characterized by a temperature scale Ti = 
(Tttv f cik f\ sin tfr\^TQ^)(Q.kBg^qQy\ where Q is the sample 
volume, cl the c-axis lattice spacing and both the Fermi ve- 
locity Vf and tfr, the angle between velocities Vk and vr+q, are 
evaluated at a hot spot. 

Low-T regime (I) When T « rTo we have 

PSn^,j^\^jTl^,,li^\ (5) 

po \3rT1T2' po 5\3rTiT2' 

In regime (I) inelastic scattering becomes comparable to im- 
purity scattering when T > ■\JrT2T1. The WF law would be 
obeyed by T = extrapolations made from regime (I) but, 
close to the quantum critical point, this regime will be hmited 
to inaccessibly low temperatures. 

Intennediate-T regime (II) For rT2 « T « T2, 

p(T) '■Ti ^, T w(T) '■T2 ^. T 

= 1 - flo— + Do— , = 1 - fl2— + 02— (6) 

po 7 1 7 1 po 7 1 1 1 

where a„,b„ are positive numbers with the n = 2 terms being 
larger. a„ are given by a„ - dx{-dfo/dx)a(x)x" 
and b„ is the same with a(x) replaced by b{x) where 
a(x) = TT"' dx'[ fo(x' - x) + no(x')]x'{x' ^ + c^)"', b(x) = 
^ ' ^x'iM^' -x) + no(x')]sgnx'[j - atan(pi)]. In regime 
(II) the inelastic scattering rate exceeds t^' once T > T^. Ex- 
trapolations to r = from regime (II) violate the WF law. 
This is because the T - intercepts of p{T) and w(T) are 
not due only to impurity scattering, they also include an in- 
elastic contribution coming from the r-linear term in Eq. 
m (WFv results from inelastic electron-electron scattering, as 
studied rigorously in the context of disordered metals. |26j]) 

In Fig. [1] 6(T) = [w(T) - p{T)]/p() is plotted and extrap- 
olated to r = from fits made above an arbitrary Tmin (sup- 
posed to be the lowest measurable T) for several values of the 
tuning parameter r. For small |r|, Tmin is in regime (II), so 
d{T) is linear in T with an r-independent slope and a non-zero 
extrapolated T - intercept that increases as r is decreased. 
Extrapolations from above r„„„ suggest T — > violation of 
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FIG. 1: Wiedemann-Franz law violation near quantum critical point. 
Left inset: Phase diagram showing critical point between a Fermi 
liquid (FL) and spin density wave (SW) with tuning parameter r. (I)- 
(III) denote dilferent temperature regimes discussed in text. Center 
inset: The thermal w(T) and electrical p(T) resistivity close to the 
critical point (r = -0.03, upper curves) and further from it (r = 0.3, 
lower curves). Main panel: The dilference 5(r) = [w(T)-p(T)]/pom 
resistivities for r = 0.3,0.1,0.05,0.03,0.01,-0.01,-0.03 from bot- 
tom to top. The supposed minimum measurable temperature is rn,i„ 
(too large to see SW or FL states for |r| << 1). The Wiedemann- 
Franz law (5(0) = would appear to be violated based on extrapola- 
tions made from above T^]„. Right inset: T = intercepts of d(T) 
obtained from extrapolations from T > Tmin, which increase as r is 
decreased. 



the WF law. The low-T recovery of the WF law, which occurs 
for r > 0, is unobservable. For r < 0, the model fails when T 
is decreased to values comparable to \r\T2 (i.e. it breaks down 
as spin-order is approached) and says nothing about T ^ Q 
WF behavior for r < 0. [Negative r values are allowed in the 
model 1 25] within regime (II).] 

The results in Fig. [T] capture much of the low-T behavior 
observed in the interlayer transport of CeCoIns. Tanatar et 
g/ ild] measured 6(T) oc T with a non-zero intercept S(T — > 
0). With decreasing field, 6(T) underwent a rigid upward 
shift and d(T 0) increased from slightly negative values 
at high fields to positive values close to He- We can make 
a semi-quantitative comparison with this data: the measured 
6{T — » 0) decreases by 0.2 as the field goes from 5.3 to 6 T. 
Using r - H/Hc - 1, taking the constant c = 1, and associat- 
ing the measured field-dependence of d(T 0) near He with 
its predicted linear dependence on r, we obtain T2/T1 ^ 10. 
Ti a: 400 mK is estimated independently from the slope 
of 6(T). So T2 ~ 4 K, which is consistent with the value 



from neutron scattering ][18[] and with the temperature below 
which WFv begins and linear T-resistivity ends in intralayer 
transport] I61]. T2/T1 = 10 and c = 1 are also used in Fig. [T] 
and the plots may be compared to those 12711 of data in Ref. 
[To. The r values in Fig. [T]correspond to H ranging from 0.2 
T below He to 3.5 T above it, roughly consistent with a range 
used to fit the T-dependence of specific heat with the same 
model. Negative r values in CeCoIns would imply that the 
r = critical field for the magnetic-paramagnetic transition 



FIG. 2: Qualitative behavior of 6{T) over a wide T-range. Main 
Panel: Solid curves at low T are plots of d(T), as shown in Fig. [T] for 
decreasing r between r = I and r = from bottom to top. The single 
solid curve at high T is the (r-independent) approximate result, and 
dashed curves are extrapolations into the crossover regime. As r is 
decreased, the peak in 6(T) narrows and shifts to lower T but the 
peak-temperature T,, remains finite at r = 0. Inset: The r-dependent 
temperature scales associated with 6(T). The peak temperature Tp, 
lies between the Fermi liquid temperature Tfl and spin-fluctuation 
temperature T2. 



He is larger than the superconducting transition field at the 
lowest measured temperatures. In the r > curves of Fig. [T] 
extrapolations made from T as low as 0.0172 ~ 40 mK would 
indicate WF violation, missing its low-T recovery. 

High-T regime (III) If we crudely extend the model to 
T > T2, then it predicts isotropic, T-linear scattering. (This 
assumes that ;t'((jj = 0) does not decrease as fast at l/T so 
x"i(^), Eq.|3] restricts lj integrals to w << kBT, which should 
be a reasonable approximation for T > T2.) We then have 
p(T) = po(l + a'JITi), and 5{T) ^ Tl/in^TTi) where a,', ^ 1 
is a constant. 

Fig. |2]is an approximate plot of 6{T) over a wide tempera- 
ture range: 6(T) increases with T before peaking at a tempera- 
ture Tp and falling off like l/T. This behavior is an alog ous to 
what is seen for electron scattering from phonons lll4ll (with 
T2 playing the role of the Debye temperature) where WFv 
is small both at low-T, where few phonons are excited, and 
high-T, where the phonon energy is small compared to ther- 
mal electron energy so scattering is elastic. As r is decreased, 

the peak in 6(T) naiTows and its position Tp shifts to lower T 

1/2 



but r,, does not tend towards zero. Rather, rp/r2 ~ (^1/^2) 
at r = 0. Similar behavior is seen in S{T) data for intralayer 
transport] 10. 16] in CeCoIns, though in the latter data the re- 
sistivity is not T-linear over the temperature range in which 
the peak occurs. 

The above assumed 2D fluctuations as suggested by NMR 
data on CeCoIns] 17]. Recent neutron data] 18] see only weak 
spin anisotropy and a dimensional crossover below IK (3D 
at low r) has been reported. ]i28i] For 3D fluctuations (as- 
suming hot-spots still exist), we write Wq^'g = {q\\ - Qy)^ + 
a^iqz - Q:)^ where Qii is in the layer and a measure of 
anisotropy, and take]]!!] {^qoY^ = 1 + c(T/T2)^'^. We find 
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P(T)Ipq = a^o[T'^^l{TiTl'^)-{rl2c)T-'l^l{TiT-'l^)] and 
w(T) = (9/5)p(r) with fl3Z3 * 1.05/(ac'/2). So resistivities 
vary as J^^^ and the field-dependent term is proportional to 
r~'^^. A 3D-2D crossover would be diflicult to observe in 
transport given the slight change in temperature power laws 
but careful analysis of the T-dependent WFv could reveal the 
dimensionality of the scatterer. 

Quantum oscillations are seen ll2ol l29ll in CeCoIns for mag- 
netic fields H as low as 3i/c, so orbital eff'ects of the field 
might be important near He- We have studied these using a 
Boltzmann-equation for interlayer transport with a field along 
k,. At T « T2, near He, hot-spot scattering is important 
since it reduces electron density in cold regions that follow 
hot spots, in the sense of cyclotron motion. One can solve 
the Boltzmann equation in cold regions, treating the density 
at the hot spot as a boundary condition obtained by solving 
the equation within the narrow range, say -(50/2 < (f> < d(p/2 
where f '(0) operates. Eqs. |5]and|6]are valid if the cyclotron 

frequency = eHvj/kf satisfies ojc > J^g^^2^^^^'^'^' "^^^^ 
is why Ti, the temperature at which the 0-averaged f"'(0) 
equals t^', is a relevant scale. The much lower T at which 
f '(0) - should not matter near He- The field helps pre- 
vent the onset of an anisotropic electron distribution ! 15h as T 
is increased. 

In summary, the field-tuned Wiedemann-Franz violation 
seen in extrapolated linear intercepts of interlayer resistivity 
in CeCoIns is explained within a Fermi liquid picture in which 
two-dimensional antiferromagnetic spin fluctuations are the 
main source of inelastic scattering. This picture naturally pro- 
duces multiple temperature scales associated with transport 
and shows how these scales vary near the critical point. 
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